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Abstract 

A decrease in current with increasing voltage, often referred to as negative differential resistance 
(NDR), has been observed in many electronic devices and can usually be understood within a 
one-electron picture. However, NDR has recently been reported in nanoscale devices with large 
single-electron charging energies which require a many-electron picture in Fock space. This paper 
presents a generic model in this transport regime leading to a simple criterion for the conditions 
required to observe NDR and shows that this model describes the recent observation of multiple 
NDR's in Spin Blockaded transport through weakly coupled-double quantum dots quite well. This 
model shows clearly how a delicate interplay of orbital energy offset, derealization and Coulomb 
interaction lead to the observed NDR under the right conditions, and also aids in obtaining a 
good match with experimentally observed features. We believe the basic model could be useful in 
understanding other experiments in this transport regime as well. 

PACS numbers: 
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I. INTRODUCTION 



The recent observation of current suppression due to "Pauli spin blockade" in weakly 
coupled-double quantum dots (DQDs) represents an important step towards the realization 



and manipulation of qubits [l|, |2|, y, 14, M • It is believed that the "Pauli spin blockade" arises 
from the occupation of a triplet state [l] that is filled, but is not emptied easily. Our model 
supports this picture and puts it in the context of what we believe is a much more generic 
model involving a "blocking" or "dark" state Q]. 

A decrease in current with increasing voltage, often referred to as negative differential 
resistance (NDR), has been observed in many electronic devices such as degenerately doped 
bulk semiconductors [sj], quantum wells 9] and even nano-structures [lol |. all of which can 
usually be understood within a one-electron picture. However, a number of recent experi- 
ments have reported NDR in nanoscale devices with large single-electron 
charging energies which may require a many-electron picture in Fock space. Although spe- 



cific models such as coupling to vibronic states 



contact couplings 



spin selection rules 1J, ll5[ asymmetric 



6|, internal charge transfers |l2j, and possibly conformational changes 
have been put forth to explain some of these observations, we are not aware of any generic 
models comparable in clarity and generality to those available for the occurrence of NDR in 
the one-electron regime. 

The objective of this paper is to present such a model for the strongly correlated transport 
regime, within the sequential tunneling approximation 15|, Il6] , where current collapse and 
NDR arise from the system being locked into a "blocking" many electron state that can only 
be filled from one contact but cannot be emptied by the other. Once it is occupied, this 
state blocks any further current flow. This concept (also invoked in a past work 7|) can be 
used to understand a general category of experiments that involve NDRs. The basic idea 
behind current collapse is a condition under which an excited state of a Coulomb Blockaded 
channel, normally inaccessible at equilibrium can, under transport conditions, be occupied. 
A drop in the current at a higher bias occurs if this excited state is a "blocking" state 
charecterized by a slow exit rate. In the first part of this paper we hence obtain a simple 
criterion for NDR in terms of the rates of filling and emptying of states. 

In the later part of this paper, we then apply our model to a specific example in 
detail, namely the NDR observed in weakly coupled-double quantum dots (DQD's) 
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Our analysis, directly addresses current collapse in terms of current magnitudes, without 
invoking non-equilibrium population analysis |l7[ fis| . Thus current magnitudes at various 
bias voltages can be solely expressed in terms of the DQD electronic structure parameters 
and electrode coupling strengths. It then becomes very transparent as to how the interdot 
orbital offset, hopping, onsite and long range correlations dramatically affect the observed 
NDR's, and leakage currents. In fact, all the non-trivial features in spin blockade I-V's [l| 
such as multiple NDR's, bias asymmetry in current, gateable current collapse and finite 
leakage currents can all be explained using the basic ideas developed here. We believe that 
our criterion obtained from a simple model could be useful for a wide class of experiments 
beyond the DQD structure discussed in this paper. 



II. CURRENT COLLAPSE MECHANISM 

Consider a generic Coulomb Blockaded system shown schematically in fig. la that is cou- 
pled weakly to contacts defined by chemical potentials /i^ and respectively. Such a 
Coulomb Blockaded system is strongly interacting and is described by its states in Fock 
space. A general condition for NDR to occur can be cast in terms of three Fock space states 
| A), | B), and | C) with energies Ea < E B < E c respectively, shown in fig. lb, that repre- 
sent three accessible states within the bias range of interest. Typically, | A), | B) could be 
the ground states of the N — 1 and N electron systems, while | C) , the first excited state of 
the N electron system. This description is equally valid for similar transitions between iV 
and N +1 electron states. Transport of electrons involves single charge removal or addition 
between states | B), \ C) and | A) that differ by an electron, via addition and removal 
transition rates Ra^b,c- Such an electron exchange is initiated when reservoir levels are 
in resonance with single electron transport channels €\ = Eb — Ea and e 2 = Ec — Ea re- 
spectively. The I-V characteristics shown schematically in fig.lc, of this three state system 
shows two plateaus with current magnitudes I p i and I p2 respectively. Current Collapse or 
NDR occurs when I jA > I p2 . Let us now derive a condition for NDR to occur in the positive 
bias I-V characteristics, of this three state system. 

These plateau currents can be evaluated in terms of bias dependent probabilities of the 
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FIG. 1: a) A generic mechanism for current collapse, rectification and leakage currents through a 
Coulomb Blockaded system can be cast in terms of b) three device Fock-space states, and transitions 
| A) <->| B) and | A) *-*\ C) between those that differ by a single electron, c) Typical I-V 
characteristics of this system comprises of two break points following a sequential access of two 
transport channels ei and ei- Insets depict the bias configurations that corresponds to various 
plateau currents. The central result of the paper, i.e., the condition for current collapse to occur 
(Ipl > Ipz), under positive or negative bias conditions, derived in Eqs. [6] and is dictated by the 
rate of de-populating the "blocking" state | C) . Such transitions are represented via dotted arrows 
in b). 



three states | A), | B) and | C), given by a set of coupled master equations 
dP A 
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24j: 



dt 
dP B 

dt 
dP c 

dt 



— {Ra-^b + Ra^c) Pa - 
—Rb^aPb + Ra^bPa 
—Rc^aPc + Ra^cPa, 



Rb->aPb + Rc^aPc 



(1) 



where Ri^j denotes the transition rate between state | /) and | J), which could be an 
addition or a removal transition between these states that differ by a single electron. These 
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rates are most generally described as 

Ra^b = 7 a M^/ Q ( ei ) 

a=L,R 

R A ^c = J2 TMl c f a {e 2 ), (2) 

a=L,R 

where j L ' R denotes the coupling strength to either contacts, M's denote coherence factors 
and fi,R(E) = f(E — ^l,r) denote Fermi-Dirac distributions in the leads. The coupling 
strengths are often expressed in terms of the tunneling Hamiltonian matrix elements Tfc Q m , 



where ka represents the k th eigen mode in a = L, R either contact [25(, and m denotes 
a quantum dot state as: r y a = 2irJ2k I T ka,m | 2 5(E — e^). Note that for the downward 
transitions the /'s are replaced by 1 — /'s. The coherence factors depend significantly on the 
structure of the many-body states involved in the transition and will be considered while 
dealing with specific examples in the following sections. We also note that Rb^c = Rc^b = 
0, since transitions between states with equal electron numbers are forbidden, in the absence 



of coupling to radiation fields {(J, or higher order tunneling 



161 ] . The steady-state current I 



equals both the left and right terminal currents (Il,Ir) given by 

q 2 L L 
Il = -^-(R a ^ b Pa - R b ^ a Pb 

+ R A _ C P A - Rc^aPc) 

<? 2 R R 

Ir = --^(Ra^bPa - R b ^aPb 

?R d t~>R 



+ R«_> C P A -R«^ A P C )), (3) 



where q is the electronic charge and Rj^j, stands for the left or right contact contribution 
to the transition | /) J). Positive bias plateau currents, shown schematically in fig.lc, 
(see appendix I), are given by: 

r Q 2 IabIba 
pl h l L AB + 7 & 

r g 2 lAB + lAC {A , 

where 

ifj = ^Mfj, (5) 
When I P 2 < I p i, current drops once the transition | A) C) is accessed. Using the above 
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expressions for I p i and I p2 , this leads to a simple condition: 

1 1 1 

"IT > —L~ + -R-> > M ( 6 ) 

7c*a Tab 7ba 

for current collapse or NDR in the I-V characteristics to occur along the positive bias. A 
similar condition for negative bias NDR is obtained by replacing L with R and vice versa: 

4 > 4 + i' ( " B> " 1) ' (7) 

The positive bias leakage current is now dictated by 7^, which corresponds to the rate 
of leakage from the blocking channel | C). For L [A > I p2 , currents rise just like a regular 
Coulomb Blockade staircase due to access of a second transport channel. Asymmetry of this 
blockade leads to current rectification observed in experiments is dictated by the ratio: 

_ I p2 (V D < 0) 

r " uy» > 0) (8) 

Thus current collapse can occur along both bias directions, provided Eqs. El [7] are simul- 
taneously satisfied. Finally, asymmetric I-V's result when Eq. [8] is satisfied, as shown in 
fig. 2b. The condition for asymmetry to occur is just based on the ratio between posi- 
tive and negative bias plateau currents, and generally depends on various couplings. Thus, 
Eq. El dand El summarize the central concept of this paper. We will show in the next 
section that they provide a clear explanation for experimental observations reported for a 
weakly-coupled DQD system. It is worth noting that these conditions can be generalized 



even if the Fock space states | A), | B) and | C) have specific degeneracies 26( provided 
the rate constants jj^' R are modified by appropriate degeneracy factors. At this point it is 
also worth noting that the degeneracies considered here are such that individual processes 
that are involved in the transitions are uncorrelated. This justifies the use of rate equa- 
tions ignoring off-diagonal effects. Inclusion of off-diagonal terms, via a density matrix rate 

not considered here. 



equation 



271 ] may be required to capture other subtle degeneracies 
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Before, we discuss coupled quantum-dots, it is useful to look into a simpler example 



discussed by various works [la . |28| and see how it maps onto our generic "blocking state" 
model (fig. lb). Consider a single quantum dot with two non-degenerate levels, one up-spin 
and one down spin, as shown in fig. 3a. The Fock space picture in fig. 3b looks just like fig. lb 
with I A) being the state with no electron | 00), | B) being the state with one up electron ||), 
I A) being the state with no electron | 00), and | C) being the state with one down-electron 



1 1 1 

— <— +— 




V D (artfLts.) 



FIG. 2: I-V characteristics: a) Forward bias I-V characteristics of the three state system with 
and set to 1 meV. As Jqa is varied, the plateau currents I p \ and I p 2 resulting from a serial 
access of transport channels e\ and €2 with bias, results in a staircase I p \ > I p 2, saturation I p \ = I p 2 
or NDR I p 2 < I p i- Forward bias NDR is achieved when Eq. [6] is satisfied, say when 7^ = 0.3meV. 
The magnitude of this blockade current I p 2 is then governed by 7^. b) Reverse bias NDR due to 
the condition (Eq. [7]), is governed by the removal rate 7^4- 



||). Note that the state with two electrons is assumed to have too high an energy to be 
accessed in the bias range of interest. It is easy to see that for positive bias that = 7^, 
7^ c = 7f, j^b = 7^, and j^c = if ■> where if 'f denotes the contact couplings with up and 
down spin reservoir channels. Ordinary contacts typically have 77 =74, so that the NDR 
criteria (Eqs. ElEJ) are never satisfied. However, if one contact, say the right is spin polarized 
7f << 7f then NDR criteria is met for positive bias though not for negative bias. 



III. PAULI BLOCKADE IN WEAKLY COUPLED DOUBLE QUANTUM DOTS 

Our idea of blocking states involves an inherent asymmetry within the transport problem. 
In our previous example (fig.3), this was achieved by explicitly postulating a spin polarized 
contact as the asymmetry mechanism. Asymmetries leading to blocking states, can also arise 
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FIG. 3: a) A simple physical realization of the above NDR condition using a single interacting 
quantum dot with spin split levels coupled to two reservoirs, one (say the left) that is unpolarized 
and the other (right) one spin polarized in the up direction, b) This situation automatically 
satisfies the NDR condition with ||) acting as the blocking state | C), as shown in the Fock space 
diagram, c) An equivalent transport energy diagram simply implies that the conducting state 
transport channel is lifted out of the bias window by an amount U due to interaction with ei 
blocked inside the device, since 7^ « 0. /iL >>| ej, ± | ensures that the blocked electron is not 
reinjected to the left contact. 

from asymmetries within the system's internal degrees, thus requiring no apriori conditions 
on bare physical contact couplings. The asymmetry effects now enter the model through 
coherence factors mentioned earlier in Eq. The rest of this paper focuses on how these 
coherence factors could lead to the NDR conditions Eqs. El [7] derived in the earlier section. 

We consider, a DQD system shown in fig. 4a, with orbital energies e± 7^ 62, and specifically 
ei > e 2 on each dot. We will now show that, for specific orbital offsets Ae =| e 2 — ei |, its 
interplay with hybridization t, and Coulomb repulsion U mp can result in a current collapse 
following the notion developed in the past section. The I-V characteristics discussed in this 
section, closely match the experimental features [l| that indicate Pauli spin Blockade. 
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FIG. 4: Many electron spectrum of a DQD system, a) Orbital energies of the two spin degenerate 
ground levels a, a and b,b on each dot, are given by e\ and €2- On-site and long range Charging 
parameters are given by Un, U22 and U12 respectively. The interdot hopping parameter is t. b) The 
many-electron spectrum comprises of 16 states as shown in Fig. 4b, labeled as | N,i). States within 
a charge subspace are just numbered and are not necessarily in any energetic order. The N = 1 
and N = 2 spectrum are relevant to the present discussion. The three-state system responsible for 
spin blockade has degeneracies and consists of two bonding | B), and one singlet | S), and three 
triplet I T) respectively, c) The first two transitions involved in the spin blockade NDR are the 
two transport channels which can be recast in the familiar tristate form shown (fig. lb). 

A. Electronic Structure of a Double Quantum dot 

In this paper, we consider a double quantum dot as shown in Fig.4a, that comprises of 
two dots whose lowest spin degenerate orbitals are labeled a, a and b, b respectively. Our 
analysis is limited to one orbital per quantum dot thus restricting our analysis to a few meV 
along either bias directions. Higher orbitals on each quantum dot can in a similar way result 
in multiple NDR's [l| along either bias directions. These higher orbital effects are accessed 
at a much higher bias than those considered here. We believe that the higher bias features 
are an extension of the general idea developed here, which will be addressed in a future 
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work. Let us now define the DQD Hamiltonian [17|, [19|, l2Cj in a localized orbital basis as: 

Hq D = ^ Cm ™ m + ( tm P C ln C P + c - c ) 
m p 

U m ^ m ^n m ^n m [ -\- — ^ ^ U m pTi m np, (9) 
p 

where the left and right dot single particle states are | a), \ a) and \ b), \ b) respectively, 
a, b denoting the opposite spin state. t mp denotes the hopping term between the two dots, 
and Urn\mi , and U mp denote the on-site and long-range Coulomb repulsion term respectively. 
The many-electron spectrum shown in fig.4b comprises of 16 states labeled as | N, i) defined 
by the number of electrons iV and the excitation number i. Six states that are relevant 
to our spin blockade discussion are labeled within the boxed sections in fig.4b, and form 
a subset of the N — 1 and N = 2 blocks. The six states consist of a doubly degenerate 
bonding | 1,0) state in the N = 1 subspace: 

| B+) = £ | a) + k | b) 

\BJ)=£\a) + K\ b), (10) 

a singlet | 2, 0) state in the N = 2 subspace: 

| S) = a (| ah)— | db)) 

+ P\aa) + 8\ bb), (11) 

and a three- fold degenerate triplet | 2, 1) states also in the N = 2 subspace given by: 

I T +1 ) = | ab) 

I T ) = -^(| ab)+ | ab)) 
I = | ab). (12) 

Here K,£,a,(3, and 5 represent the wavefunction coefficients for the various basis Fock- 
states. We shall see in the following subsection that these wavefunction coefficients depend 
on the physical parameter set £1,62,^^11,^22,^12 of the DQD, and hence appear in the 
calculation of various coherence factors (Eq. [2]) involved in the transition rates. It is worth 
noting that the anti-bonding level | 1,1) of the N — 1 spectrum is not involved in the 
Pauli-Blockade NDR. This is because in the weakly coupled DQD that we consider here, 
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the bonding-antibonding gap is directly proportional to the orbital energy offset Ae, which 
is of the order of a few milli- volts. The energy scale relevant to our NDR discussion is the 
singlet-triplet splitting which is of a couple of orders of magnitude smaller given that the 
DQD system is weakly coupled. The anti-bonding state then becomes relevant beyond the 
bias range considered here where the spin-blockade is lifted. 

Coherence Factors: Although six states, (fig.4b) are involved in the spin blockade I-V 
characteristics within the bias range of interest, they can be cast in the familiar tri-state 
form (fig. lb) as shown in fig. 5, by incorporating the appropriate degeneracy factors j3] for 
the transition rates. The degenerate states may just be lumped as one state under such 
specific conditions 2fJ. The coherence factors are evaluated between two states that differ 
by an electron, whose structure and symmetry properties depend on the physical parameter 
set of the DQD and how the individual dots couple to the two electrodes. The general 
definition is as follows: 

^ = l a \(N,i\cl\N-l,j)\ 2 , (13) 

where is the creation/annihilation operators for an electronic state on the end dot coupled 
with the corresponding electrode, and 7° is the bare left or right electrode coupling factor 
defined in section I. In our case, the coherence factors are evaluated between the a) | 1,0) 
bonding and | 2,0) singlet: 

Too = 1 L I (S I I B.(B+)) | 2 
= 7 L (^ + "k) 2 

loo = 1 R I (S I I B.(B + )) | 2 

= 7 >£ + k<5) 2 , (14) 

b) I 1,0) bonding and | 2, 1) triplet states: 

7i L o = 7 L (I (To I c{ I B.) | 2 

+ I (T I 4 I B+) | 2 ) 

= I (T±i I < (B) I B T ) | 2 = 7 V 
7fo = 7 B (I (To I 4 I B.) | 2 

+ I (To I 4 I B +) I') 

= l(T ±1 |ct (6) | 5T )| 2 =7^ 2 . (15) 
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Recall that the transport feature central to the Pauli Blockade experiments is the occurrence 
of multiple NDR's l| along both bias directions, and this condition must be recast into the 
model developed in the first section. Therefore, the basic criterion for NDR to occur along 
the positive bias (Vp > 0) direction is: 



1 1 1 

27fo 7o L o 2 7oo 



>— + ( 16 ) 



and along the negative bias direction (Vz> < 0) is: 

1 > 1 + 1 (17) 
^7ib ^loo 7oo 

Notice that NDR conditions in the above equations also include degeneracy factors as de- 
picted in fig. 5. A direct substitution from Eqs [T^l and [15] with j L = ^ R yields 

111 

2? > + aKf + 2{ai + k5) 2 (18) 

along the positive bias direction and 

1 1 1 

3^ > 2(/?e + olkY + « + k5) 2 ' ( } 

along the negative bias direction. To apply these relations we need to consider the wave- 
function coefficients £,K,a,(3, and 5 in more detail. 



B. Spin Blockade Transport 

Consider first the wavefunction coefficients k and £ for the N = 1 bonding state | B±). 
It is easily seen that k ~ £ for Ae << t and k >> £ when Ae >> t, given that e± > e 2 . The 
two electron triplet (see fig. 5) is relatively un-affected by electronic structure variations, 
but the singlet state described by the wavefunction coefficients a, (3 and 5 requires some 
more discussion. 

The presence of a single electron (fig. 6) in one dot influences the addition of another on 
either of them as a result of on-site or long range Coulomb interaction. For t << Ae, U 
relative magnitudes of a, (3 and 5 now depend on off-sets between the total energies of the 
three configurations shown in fig. 6 which we shall consider now. 

Transport Scenarios: Consider three transport scenarios depicted in fig. 7, which are based 
on the relative offsets between orbital energies Ae = e x — e 2 , and relative offsets between 
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N=2 




FIG. 5: Internal structure of the Three state system: This three-state system responsible for spin 
blockade comprises of one-electron Bonding | B), and two electron singlet | S), and triplet | T) 
respectively. Electronic configuration of the three states is also shown along side. The singlet, 
triplet and bonding states form the trio of Fock space states involved in the SB I-V characteristics, 
with the triplet state being the blocking state. 

the two electron wavefunction configurations (fig. 6) AEi — t\ + Uu — €2 — c7 12 , and AE 2 = 
£2 + U22 — ei — Uu- In all our simulations we have used parameters consistent with a 



weakly-coupled DQD structure in jl], with t = 0.2 — 0.3meV, Uu = U22 = 4meV, and 
j L = j R = O.OlmeV, while other parameters are varied to the three cases depicted in Fig. 7. 

Case I: Coulomb Staircase: In first scenario shown in fig. 7a, Ae = 0, Uu = U22, which 
implies k = £ = ^= and (3 = 5. A straightforward evaluation of coherence factors from 

HI [T7| is never satisfied, leading 



30j. In this case, both left and 



Eqs [15] implies that the criterion for NDR given in Eq. 
to regular CB staircase I-V characteristics shown in fig. 8a [c 
right contacts act are equally efficient in electron-addition and removal processes and hence 
no such blocking states can be formed. Thus the NDR phenomenon is absent. 

Case II: Forward bias Current Collapse: In this case 2e2 + U22 ~ £1 + £2 + U12, which 
implies that the two electron state | bb) is almost degenerate with (-^(l ab)— \ ab)) implying 
that both electrons, whenever in a singlet state can either reside inside dot 2 alone or one 
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FIG. 6: Understanding the internal structure of the 2 electron singlet: Here, both on-dot U mm 
and long range Coulomb repulsion U mp affect the wavefunction, by raising the onsite energy of 
the second electron depending on where the first electron is localized. The singlet state is a 
superposition between states with one electron in each dot (| ab), \ ab), \ ab),\ ab)), and the ones 
with both on either dots (| aa) and | 66)). A remarkable consequence of the interplay of electron 
correlations, hopping, and orbital energy offset is the fact that the singlet state is a mixture of 



states localized on the same dot with states present on both \2W . The wavefunction coefficients a, 
f3, and 5 signify the contribution of each configuration and can be varied by tuning the DQD's 
internal structure parameters. 

in both dots. The condition e\ > 62 makes it energetically unfavorable to access | aa) 
resulting in a ~ 5 >> (3. Under these conditions, note that our forward bias NDR condition 
(Eq. [E2 dHD results in 



2 > 7-^2 + ^2 ( 2 °) 



1 1 1 

2e > ( Ka ) 2 + 

which is trivially satisfied in our case e\ > €2- The I-V characteristics now show a prominent 
NDR as shown in fig. 8b, noted in the experimental trace shown in fig.8d. The parameter £ 
dictates how effective the right contact (also see Eq. [T5j) is, in the electron removal process 
between the triplet and bonding states. Recall that £ << k when e\ > 62, implying that the 
right contact has a very in-efficient triplet removal, thus blocking the triplet within the DQD 
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FIG. 7: Transport scenarios through a DQD system, a) Zero orbital offset Ae = 0: When XJ\\ = 
U22 > U12, NDR condition is not satisfied under any bias conditions, b) Finite orbital offset Ae > 
with resonance e\ + U12 = £2 + ^22 : Condition for positive bias NDR is satisfied, c) Finite orbital 
offset Ae > and off-resonance e\ + U\% < €2 + ^22 : the NDR condition is satisfied for both bias 
directions. 

system. This causes a current collapse once the population of triplet state is energetically 
feasible. The Forward bias I-V shown in fig. 8b shows a strong NDR with current collapse 
leading to a small leakage current, once the bias permits the access of transport channel 
eio in the bias window. In the Fock space picture (fig.4c,5), current at lower bias due to 
I B) S) drops once the transition | B) T) occurs since | T) is a blocking state. 
Note that, when transitions \ B) *-*\ S) and | B) <r+\ T) are simultaneously accessed, the 
system is put in a blocking state once threshold is reached, thereby no NDR occurs. This 
is shown dotted in fig. 8b and occurs experimentally [2], upon application of gate potential, 
thereby changing the relative position of threshold voltage. We also find that, a finite but 
small leakage current in the order of 5 — 10 pA occurs due to a finite 7^ oc £ 2 corresponding 
to the leakage of triplet into a bonding state. Importantly this leakage current reduces by 
further detuning the dots, making the bonding state more localized on the second dot, such 
that 7f oc £ 2 — > 0. Since V = 7^7^ > 1 we see a clear asymmetry between the two bias 
directions in the I-V characteristics. 
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FIG. 8: I-V Characteristics corresponding to various transport scenarios: a) Regular Coulomb 
Blockade plateaus in the absence of orbital off-set, since NDR condition not satisfied, b) Forward 
bias NDR occurs when there is finite orbital off-set along with resonance condition corresponding 
to fig. 7b. The NDR condition derived in Eq. [6] is physically achieved in the DQD system by 
orbital offset such that k » £. Furthermore, the resonance condition U22 — U\\ = e% — e 2 implies 

a = 5, thus resulting in -4j- >> -L- + -4-. Also shown dotted is the fact that NDR need not occur 

7io 7 00 Too 

when the blocking state is simultaneously accessed with the conducting state. This can can occur 
by varying gate voltage Q]. c) Reverse bias NDR's can also occur in the presence of finite orbital 
off-set with off-resonance discussed in fig. 7c. d) and e) Experimental trnces Q] corresponding to 
scenarios b) and c). 



Case III: Multiple Current Collapse: Seldom do parameters exactly match the condition 
2e 2 + U22 ~ ei + e 2 + U12, experimentally thus making it possible for reverse bias NDR's 
too. This corresponds to the scenario depicted in fig. 7c. The important consequence of 
Eqs. E], [7] was that NDR along a bias direction only depends on how the blocking state is 
coupled to the emptying reservoir. Now the negative bias (Vd < 0) NDR condition simply 
is -^r- > -jt + ttt, which is also satisfied when e x + Lq 2 < e 2 + U 2 2- Under these conditions, 

d 7io Too z ~too 

shown in fig. 7c, 8c a mild NDR is noted even in the reverse direction, in good match with 
experimental (fig.8d,e) observations 

The experimental I-V's also show a lifting of spin blockade followed by successive NDR's 
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at higher positive bias directions. The lifting of NDR's can be explained simply due to access 
of transport channels between two electron higher excitations and one electron antibonding 
states. A consistent match with experimental features, will however involve the inclusion of 
more orbitals within each dot, such as p-type orbitals, which we shall address in a future 
work. The fundamentals of NDR's, collapse and rectification still remain within our general 
framework derived in this paper. 

IV. CONCLUSION 

In this paper, we have developed a general model for multiple current collapse (NDRs) 
to occur in the I-V characteristics of strongly interacting systems using the idea of blocking 
transport channels. With the aid of this model, we have provided an interpretation of non- 
linear transport through weakly coupled quantum dots in the Pauli spin blockade regime. 
Our model captures subtle experimental observations in this regime, that includes multiple 
current collapses which are gateable, leakage currents, and rectification, all of which are 
consistent with the experimental double quantum dot setup. We believe that this model 
can be exploited in the understanding of other novel blockade mechanisms that can arise 
from different internal degrees inside a Coulomb Blockaded system. 

Acknowledgments: It is a pleasure to acknowledge Prof. Avik Ghosh for stimulating 
discussions. This project was supported by DARPA-AFOSR. 

V. APPENDIX I : DERIVATION OF PLATEAU CURRENTS 

The plateau currents I p \ and I P 2 shown in Fig. lc are the saturation currents due to 
transitions | A) B), and | A) C) whose one-electron energies are ei = Eb — Ea 
and e 2 = E c — Ea respectively. For plateau currents, (say) under positive bias > Hr 
shown in fig. lb), additive transitions rates are given by Ra^(b,c) — ^a-^(bc) anc ^ remova l 
transitions are given by R(b,c)-*a = Rfsc)^A- Furthermore, in order to evaluate I p i, given 
a sequential access of the two transitions, only states | A) and | B) appear in Eq. [TJ and El 
A steady state solution of this master equation Eq. [TJ with P c = along with Pa + Pb = 1 
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gives: 

Pa 



Rb^a 
Ra-^b + Rb^a 



PB = R ■ (21) 

Using Eq. [3l and the fact that f L (ei) = 1 and f R {ei) = for positive bias, plateau current 
I p i is expressed as: 

Ipi = ~r {Ra~>bPa) 
= (R^aPb) 

= q 2 Ra^ b Rb-,a ^ 22 ^ 
h Ra~>b + Rb^a 

At a higher bias, the transport channel e2 = Eq — Ea, due to transition | A) C) is 
accessed and thus we have to evaluate the stationary solution of the three state master 
equation Eq. [TJ The steady state solution of Eq. [TJ along with Pa + Pb + Pc — 1 yields the 
non-equilibrium probabilities of the three states as 



Rb-*a Rc-,a 

Ra->b 1 

~ i2 B _ A l + W + 

Rb~>a Rc->a 

Pc = rT-Ii+^+P^ (23) 

tiB^A K-C^A 

Similarly as in Eq. [221 with the fact that f L {t2) — 1 and f R {^2) — 0, the second plateau 
current l v i is expressed as: 

q 2 

I P 2 = ~r (R a ^bPa + Ra^c p a) 
= ^ (R^ a Pb + R^aPc) 



q 2 Ra-,b + r a-,c 

H 1 I R A^B I ^A—C 

B — >A n C^A 



(24) 



The expressions derived above are cast in terms of rates in general. One has to observe that 
each plateau current corresponds to a situation when only one contact contributes to the 
addition of an electron while the other to the removal. This implies J L (e) = 1 and f R (e) = 
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for positive bias (/i^ > /j,r) and vice- versa. Hence the rates, using Eq. [2] are now given by: 



R L A _ 




= 1 L M L AB 


= lAB 


R A- 




= 1 L M L AC 


= 7 AC 


r>R 
ti B _ 


->A 


= 1 R M% B 


= lBA 


r>R 
ti B _ 




= J R M% B 


= lCA 



(25) 

Under these conditions it is easy to see that I p i and I P 2 follow the expressions in Eq. HI and 
that conditions for NDR's under positive and negative bias conditions (Eqs. MU§ are easily 
obtained by replacing 7^'s by 7^'s. 
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